Semiclassical theory of the photogalvanic effect in non-centrosymmetric systems 
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We develop a semiclassical theory of nonlinear transport and the photogalvanic effect in non- 
centrosymmetric media. We show that terms in semiclassical kinetic equations for electron motion 
which are associated with the Berry curvature and side jumps give rise to a dc current quadratic 
in the amplitude of the ac electric field. We demonstrate that the circular photogalvanic effect is 
governed by these terms in contrast to the linear photogalvanic effect and nonlinear I-V characteris- 
tics which are governed mainly by the skew scattering mechanism. In addition, the Berry curvature 
contribution to the magnetic-field induced photogalvanic effect is calculated. 

PACS numbers: 



I. INTRODUCTION AND 
PHENOMENOLOGICAL DESCRIPTION 

Recently, considerable progress was made in the semi- 
classical theory of the anomalous Hall effect (see the 
review article, Ref. [l|). The renewed theory based on 
the generalized semiclassical Boltzmann equation re- 
moves the previous controversies surrounded this effect. 
This theory takes into account the following mechanisms 
of the anomalous Hall effect: (i) the skew scattering 
contribution which appears on going beyond the Born 
approximation^i^, (ii) the contribution to the electron ve- 
locity due to the Berry curvature in the electron Bloch 
wavefunction^i^!^ and the contributions of side-jumps, 
i.e., coordinate shifts at the scattering eventSfii^ (iii) to 
the electron velocity and (iv) to a change of the electron 
energy upon a scattering in the presence of an exter- 
nal electric field. In the light of this advance in linear 
transport physics, the semiclassical approach to nonlin- 
ear transport should be revised as the next natural step. 
In the present paper we develop a semiclassical theory 
of a dc electric-current response quadratic in an external 
dc or ac electric field. The former response is a non- 
linear contribution to the I-V characteristics while the 
latter can be interpreted as a low-frequency photogal- 
vanic effect. The microscopic theories of photogalvanic 
effects induced by light waves of sufficiently high frequen- 
cies lying in the terahertz, infrared or visible regions are 
based on quantum-mechanical calculations. They make 
allowance for both the skew and side-jump mechanisms 
(usually called respectively the ballistic and the shift con- 
tributions) and have no need for the introduction of the 
Berry curvature, for details see the books Refs. mini. 
On the other hand, the existing theories of static cur- 
rents quadratic in the electric field regard the skew con- 
tributions only, see Refs. fl^fTsl . The shift contribution to 
the I-V characteristics has been considered in the regime 
of streaming achieved in strong electric fields As for 
the Berry-curvature mechanism of nonlinear transport, 
to the best of our knowledge, it has not been analyzed 
yet. 

Phcnomenologically, the electric current density j 



which is quadratic in the electric field is described by 

jx = ^XtiuEi_iE* + Txfj.urjE^E*qj^ , (1) 

where A and T are material tensors, E and q are the 
electric-field amplitude and the wave vector of the ac elec- 
tromagnetic wave, Xyfijiy and rj are the Cartesian coor- 
dinates. In the particular limit of a static field (the wave 
frequency w — > 0) the amplitude E reduces to the dc elec- 
tric field, and q ^ 0. The second term in the right-hand 
side of Eq. (P) describes the photon drag effect. It was 
first observed and described as the high-frequency Hall 
effect by Barlow In the quantum-mechanical descrip- 
tion of optical transitions, the drag effect arises due to 
momentum transfer from photons to charge carriers pi^jlS 
This current can be induced in both centro- and non- 
centrosymmetric systems. The semiclassical theory of 
this effect is well developed, see, e.g., Refs. I20II2TI and 
the review Ref. [l^, and here we focus completely on the 
first term in the right-hand side of Eq. It can be 
conveniently decomposed as 

n - X,,AE^.E:}, + 7Am i {E X E*)^ , (2) 

where {Ef^E*}s is the symmetrized product {Ei^^El + 
E*Ei,)/2. The two contributions in Eq. ^ describe, re- 
spectively, the linear and circular photogalvanic effects, 
or LPGE and CPGE,2iiOiii At w ^ 0, the vector E is 
certainly real and the CPGE vanishes while the LPGE 
reduces to the quadratic conductivity. Since the current 
density is a polar vector and, therefore, odd under the 
operation of space inversion, while the square of the elec- 
tric field is even, the CPGE and LPGE are allowed only 
in materials that lack inversion symmetry, respectively, 
in gyrotropic media and piezoelectrics. 

We assume that, in addition to the dc or ac electric 
field, a static magnetic field JB is applied to the sample, 
so that the coefficients Xx^^ a-nd 7a^ are B-dependent. 
In the present work we keep only zero- and first-order 
terms in B, in which case one has 

7A^(B) = 7i"'+/9?^.i?.- (3) 
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Since the homogeneous magnetic field is invariant under 
space inversion the tensors f3^ and are also nonzero 
only in non-centrosymmetric media. 



part. Moreover, in the presence of a magnetic field the 
asymmetric scattering is allowed even in the Born ap- 
proximation.— 



II. GENERALIZED BOLTZMANN EQUATION 



Zero magnetic field 



Since the photon drag current is neglected we assume 
the ac electric field E{t) to be homogeneous and ignore 
effects related to the magnetic field of the electromagnetic 
wave. This means that the electron distribution function 
fp {t) is independent of the spatial coordinate r. Since 
in this paper we ignore the spin-orbit interaction this 
function is also spin- independent. The kinetic equation 
for fp (t) has the standard form 



dfp 
dt 



(4) 



where p is the time derivative of the momentum of the 
center of mass of the electron wavepacket, and Ip{f} 
is the collision term. The latter consists of two terms, 

Ip^^ {f} and /p™^^{/}, representing the elastic scatter- 
ing by static defects or disorder and inelastic scattering 
by phonons. For simplicity we assume the typical elastic- 
scattering time T to be much shorter than that for inelas- 



tic scattering, Tinei- In this case Ip^^ {f} and Ip^'^^' {f} 
describe, respectively, the electron momentum and en- 
ergy relaxation. Taking into account the Berry curvature, 
side-jumps under electron scattering, and skew scatter- 
ing, the equations for p and Ip{f} are modified as com- 
pared to the classical Boltzmann equation. Following 

Refs. IT^flsl (see also Ref. [l|) we expand the term /p°'' {/} 
into powers of the scattering potential V{r) and write 
it as a sum of the scattering rate found in the lowest, 
second-order, approximation 



r(incl)j 



E 



fp), 



(5) 



and the antisymmetric part of the scattering rate (the 
skew contribution) 

p' 

calculated beyond the Born approximation (third- and 



(a) 



(a) 



p',-p- Here Vp>p 



higher-order corrections): W^. 

is the transition matrix element, Sp is the electron energy, 
and the angle brackets indicate averaging over the ran- 
dom scattering potential. One can see that the (5-function 
in Eq. ([6]) contains a change of energy of the scattered 
electron under the side-jump rp>„ in the presence of an 
external electric field E, see Ref. Note that the inelas- 
tic collision term /p"°'^{/} also contains an asymmetric 



In the absence of an external magnetic field one has 
p = eE{t), where e is the electron charge. The electric 
current density reads 



- y^Yl ^pfp 



(7) 



Here Vd is the macroscopic volume of the d-dimensional 
sample, d = 3,2, 1 for a bulk sample, quantum well and 
quantum wire, respectively, and the factor 2 takes into 
account the electron spin degeneracy. The velocity of the 
electron wavepacket is given by 



' p 

5r. 



Up 



px Fp + 5rr. 



(8) 



p ^ ^'^p'p ^p'p ' 



where Vp = dsp/dp is the electron conventional velocity. 
The first two terms in Eq. ([8]) were derived in Refs. 
while the third term was derived in Refs. (see also 
Ref. [l] for a review). The Berry curvature Fp is defined 
by^iii^ 




^) (9) 



d / 

Fp^ — xn{p), n{p) = i{uk 



with k — p/h, and Uk{r) being the periodic amplitude of 
the Bloch wavefunction. The spatial shift of an electron 
under the scattering p ^ p' has the form^*^ 



Im ( V;,^Dp,pVp,p 



O(p') - ^l{p) , (10) 



where 



Dp'P = M + 



B. Nonzero static magnetic field 

The equation for p is modified in a magnetic field B 

(11) 



p = e 



E{t) + -^xB 
c 



where Vp is defined by Eq. ([5]). 

Taking into account the Berry-curvature-induced 
renormalization of phase-space volume^ the equations 
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for the electron and electric current densities. Nd and j, 
are also modified as follows 

Due to the same reason, in Eq. ([5]) one should perform 
the replacementi^ 

p' p' 

I 



rp = 
P = 



III. NONLINEAR TRANSPORT IN THE 
ABSENCE OF MAGNETIC FIELD 

The electron distribution function fp is expanded in 
powers of the electric-field amplitude E defined by 

E{t) = Ec-''^^ + E*e''^* . (16) 

For the classical Boltzmann equation in the relaxation- 
time and effective-temperature approximations one has 

= e--*4L^ + c.c. , = -erUE ■ Vp)f[, (17) 

for the linear correction, and 

/i,'' = foisp, 8) - Msp) + f^'\p) , (18) 
fi'\p) = 2e\Rc{T^)f[;{vp^Vp,~^^S^.){E^E:}s , 

for the time-independent second-order correction. Here 
the momentum relaxation time r is defined by 

^cu)|^|^_4_l^ ^ (19) 

where the bar means averaging over the directions of p, 
hereafter for simplicity we neglect the dependence of r 
on the electron energy Ep, /o(ep,6) is the Fermi-Dirac 
distribution function at the effective temperature dif- 
ferent from the bath temperature T due to the heating 
of the electron gas, /o(ep) = foiSp^T), a prime means 
derivative over Spi /p = dfo{ep)/dep, and 

T 

= : • 

1 — lUJT 

In the following subsections we will take into account 
skew-scattering. Berry-curvature, and side-jump effects. 



Below we take into account contributions to the current 
d?]) up to the first order in Fp and Tp/p. This allows one 
to approximate the (5-function in Eq. ([6]) by 

S[ep' -Sp- eE{t) ■ rp,p] (14) 

d 

= S{ep' - Ep) + eE{t) ■ Tp'p Q^^i^P' - £p) ■ 

Expanding Vp and p up to the first order in Fp we can 
present the solution of the coupled equations ^ and pT|) 
in the following form 



(15) 



find corrections to the electron distribution function 
given by the sum of (fT7|) and (fT5|) . substitute these 
corrections into Eq. ([7]) and find the electric current 
proportional to bilinear products E^E* in accordance 
with the phcnomcnological equation It should be 

mentioned that, for an ac electric field, the current ([2]) or 
(O is defined as the time average. Taking into account 
the definition (fTB|) of the field amplitude the tensor com- 
ponents X\ij,u {'^) in Eq. (m found in the limit lo — > +0 
and those for the static electric field differ by a factor of 2. 



A. Skew-scattering contribution 

In this subsection we set all vectors Vpip and Fp to 
zero so that the nonlinear current appears only due to 
the asymmetrical parts of elastic and inelastic scattering 

rates, respectively, Ip^'°^ {f} and /p"''''°-'{/} and one has 

^' = ^ E [if^"' {f^' ] + i^'^-Hfoisp, e)}" . 
p 

(20) 

The second term contributes to the current only in sys- 
tems possessing a polar axis. 

Equation (|20p demonstrates that the skew-scattering 
processes give rise only to LPGE. The current propor- 
tional to the degree of circular polarization of the electric 
field does not appear due to skew scattering. 

B. Berry-curvature related current 

Let us now ignore both the asymmetry of collision rate 
and side-jumps at the scattering and analyze the effect of 



v„ 1 + -B-F^ 



Eit)(l 



-B Fr 



eE{t) X Fp - -{vp- Fp) B + 5rp , 
^!lL±^xB-l(Eit)-B)F, 



p 
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the Berry curvature on the nonlinear current. Note that 
in centrosymmetric systems Fp ~ F^p, while in systems 
with time reversal symmetry Fp = — i^-_p.— Thus if the 
system is both time- and centrosymmetric the Berry cur- 
vature vanishes, Fp = 0. 

The nonlinear current is obtained by taking into ac- 
count linear-in-£J corrections to the velocity Vp and to 
the distribution function, see Eqs. p^]) and P7)) . and 
time averaging the sum 



The result written in terms of the tensors X^^^]^ ^-nd 7^^*^^ 
reads 



(0) 



(0) 



2e 



^Im (t^)^/oCa^ , 
p 



(21) 



where eA^r; is the antisymmetric unit third-rank tensor 
and 



CA^(ep) — Vpx Fp^^. 

While deriving the equation for 7^'^^ we took into account 
that the sum C^^ri makes no contribution to the cur- 
rent because Vp • [Vp x r2(p)] = 0. 

In the particular case of d-dimensional degenerate elec- 
tron gas with a parabolic dispersion Sp = /2m* {m* 
is the electron effective mass) where /cb^ ^ Ep and 
/q Ri ~S{ep — Ep) with Ep being the Fermi energy, 
Eqs. (PT|) are reduced to 

(0) em*d \Ti I \ 



2Ef 

(0) em*d 
7a^ = -^^Cx^,llTl{a^) , 

where a^^ is the linear Drude conductivity, 



(22) 



and the value of Cxfj,{e) is taken at the Fermi energy. 

Symmetry analysis allows one to find linearly indepen- 
dent coefficients in the expansion of the Berry curvature 
Fp in powers of the momentum p. To illustrate, below 
we give examples of such analysis for the bulk crystal 
symmetries T^ (zinc-blende lattice) and Cgt, (wurtzite) 
and the point symmetry of a quantum-well structure 
grown along the low-symmetry axis [hhl] different from 
[001], [111] or [110]. 

In the Td point group the components of the polar 
vector Vp (X p and axial vector Fp transform accord- 
ing noncquivalcnt irreducible representations F2 and Fi. 
Since the direct product of the representations F2 and 



Fi does not contain the identity representation Ai the 
average of the product VpxFp ^ over the directions of p 
vanishes, and the Berry-curvature related current is for- 
bidden for the Td symmetry. The zero value of 7^*^^ agrees 
with the fact that the Td crystal class is not gyrotropic 
and forbids the circular PGE. 

In crystals of the Cgi, symmetry with z being the po- 
lar axis the components Vpx,Vpy and —Fp^y, Fp^x trans- 
form according to the representation E2 and the nonzero 
averages Cx^ are C^y = —Cyx- It follows then that 
nonzero components of the photogalvanic tensors are 

Xxlx = XxL = Xyzy = Xyyz and 7^°^ = -7^°' , and they 
can be expressed via Cxy{sp) by Eqs. ([2T|) . Expanding 
Fp to first order in p we have 



Ac X p . 



(23) 



where c is a unit vector along the polar axis z and A is 
a constant. Note that a similar equation follows for the 
continuous group Ctxju- From Eq. (|23p we obtain 

and using Eqs. (|2ip we get for arbitrary degeneracy of 
the electron gas 



X. 



(0) 



A em* Re (o-^), 7^' 



(0) 



-A em* Im((T^) . (24) 



In the band-structure model of a wurtzite-typc semi- 
conductor including the conduction band Fic, and the 
valence bands r6«,ri-u we obtain the following estima- 
tion for the constant in Eq. 



A 



EgA,{Eg + A,) 



(25) 



Here 



Pi_ = —{S\px\X) = —{S\py\Y) , P|| = —{S\p:,\Z) , 
mo ma mo 

Q^^{Z\px\X) = —{Z\py\Y), 
mo mo 

mo is the free electron mass, Px ~ —ihd/dx, Eg is the 
fundamental energy gap, Ac is the crystal splitting, and 
S, {X, Y) and Z are the F-point Bloch functions in the 
Fic, Fei; and Fit, bands, respectively. Note that the polar 
crystal symmetry is fixed by a nonzero matrix element Q. 

In a two-dimensional system of the Cg symmetry with 
the interface plane normal to z and the mirror reflection 
plane (yz) a single nonzero component of F is Fp z- The 
first-order term is given by 

Fp,, ^AsPx, (26) 

where A^ is a constant. As a result one has 

Cxv{£p) = AgEpSxxS^z 

and one can readily use Eqs. ([2T|) to relate the photo- 
galvanic tensors to the real and imaginary parts of the 
two-dimensional conductivity cr^. 
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C. Shift currents 

111 the semiclassical theory, the shift current consists of 
two contributions, the first one is due to the scattering- 
induced correction to the electron velocity, see Eq. ([5]), 
and has the form 



(27) 



while the other is related to the field-induced correction 
to the collision term and can be presented by the time 
average of 



reverse under the interchange of p' and p, Vp pi = —rp/p. 
For bulk T^- and Cgn-symmetry crystals the second-order 
terms are 



rp'p = viPzP ~ p',p') 

^p'p-X miPyPz -p'yPz) 



(T.), 



(31) 



and two other components arc obtained by cyclic permu- 
tation. For simplicity, mixed terms of type PzP'\ — Pxp'z 
allowed by the Cei; point group are neglected. 

For a quantum well of the Cg symmetry, the elemen- 
tary spatial shift has the form 



m [fx - p'x) + m (pI - p'yj y 

- m {PxPy - p'xP'y) * 1 (32) 



where 



5U = -2er^[/(i)(t) - f!^}{t)] [E{t) ■ rp,p 



^{\Vp,p\')^5{e,-e,.) 



(28) 



Taking into account the identity 

dS{ep) _ dS{ep) 



dSr, 



dp 



the second contribution may be reduced into 



47re^T 



5{ep - £p') 



dp 



■(sc) 



{fi>'j-.rp'L){\Vp'pr}iE*-rp'p) 



c.c. . 



We ignore the shift-induced corrections to antisymmet- 
ric part of the scattering rate W^f^ because the combined 
effect of the Berry curvature and side-jumps is negligible. 

In order to obtain simple analytical equations for the 
photogalvanic tensors, we simplify the model and assume 
that the scattering matrix element Vp'p is the sum of a 
main term Vq independent of the electron momenta, and 
an additional term SVp'p dependent on the momenta, 
governing the spatial shift pop but making no influence 
on the relaxation time r. In this model the current j^*"^^ 
can be written as 



where x and y arc the unit vectors along the Ox and 
Oy axes, respectively. By using Eqs. (|?7)) . (|5D|) and ([5^ 
we obtain the following nonzero tensor components for a 
heterostructure of the Cg symmetry 



X^l(vl) 



(sc) 



-X^^,(vl) = 2(771- 772) em*Re(aJ , 
773em*Re(cr^) , 

2xi°Usc) = 2(771 + 772) em*Rc{a^) , 
(??i +^2)e7n*Im(cr^^) . (33) 

Since the tensor Xxli/ symmetrical with respect to in- 
terchange of the second and third indices the components 



^yyy 
v(0) 



Yx'Jisc) 



(29) xi°x(vl) and x^xyxisc) are also nonzero 



(0) 



IV. NONLINEAR TRANSPORT IN THE 
PRESENCE OF MAGNETIC FIELD 

Here we outline the derivation of the Berry-curvature 
related photocurrent linear in a static magnetic field B. 
In the presence of a magnetic field the kinetic equation 
has the form 



dfp^^,dfp 
dt dp 



1 



-B 



Fp')W^^l{fp,~fp), (34) 



where p is given by Eq. Ijl5p . It is instructive to introduce 
the auxiliary function 



,-(sc) 



2e^ 



dp 



E 



pp' 

9/0 (£p') 



V dp 



dp' 



{E* 



I p'pj 



(30) 



c.c. 



The further simplifications become possible if we expand 
the shift (fTU|) in powers of p',p and retain quadratic 
terms. The latter are the lowest nonvanishing terms tak- 
ing into account the following properties of the polar vec- 
tor Tp/pi spatial shifts under electron scattering are sym- 



metrical functions of the momenta, T"_p'._p 



fp'pj 



and 



1~-B-Fp fj 



Then in the linear-in-B approximation we have the fol- 
lowing kinetic equation for this function: 



d(j)p 
dt 



^E.^[iB-Fp)<Pp] 
c dp 



(35) 



xB -{E-B)Fp 



-(B-Fp) 



d(j)p 
dp 
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Let us introduce the notations for corrections 4>^^ ^ ^ 
to the auxihary function with to = 0, 1, 2; n = 0, 1, e.g., 
(I}^p\ 4^^^ etc. From Eq. ^ we find that the 

following expression for the electric current quadratic-in- 
E and linear-in-B 



P 

-e{E X ^) 



(36) 



-Fp X {vp X B) 



As an example we further consider a 2D system of 
symmetry and take Fp in the form Fp \\ z, Fp ^ — A^-p^^ 
see Eq. (P^ . Solving the kinetic equation ([55]) one can 
obtain the linear- and quadratic-in-£J corrections at _B = 
given by Eqs. pT]) and and linear- in-B corrections 
in the form 

0(f) = --B.^p,./o(£p) , 



r^Vp ■ {B X £;)c 



(37) 
(38) 



and, finally, 



-TT..,B,E* 



dp 



[Fp^,XE-Vp)Q+c.c. (39) 



In the above expressions, we neglected terms which do 
not contribute to the current. For the particular case of 
Cs symmetry nonzero components of the magnetopho- 
togalvanic tensors /^^^^^ and /3j^j, defined according to 



Eq. ^ become 

h'xxxz 



4/3o 



0^ 



1 + {uJTf 

4/3o 



yyz 



yyxy 



r^xy 



4/3o 



[1 + (c^t)2]2 
2 + (lut)^ 



[1 + (c^t)2]2 



L l-(c.r)2 



c 

yz2 



yxzy -^^u _^ ^^^^2]2 

c^r(l - (a;T)2 



r'xxy 



4^2 , 



[1 + (t^T)2]2 

-- -2/3oi 



1 + {lotY 



(40a) 

(40b) 

(40c) 

(40d) 

(40e) 
(40f) 



where (3q — AsN2De /m*c. 

It follows from the above equations that in the static 
limit Lu = the magnetic-field induced corrections are 
given by 



dj, = 2A,em*aH{El~E^), 

at B II z, and 
Sjx = Asem*aHEyEz, Sjy 

at B II y. Here ajj is the Hall conductivity equal to 
'-^cT'^o J where tUc is the cyclotron frequency and ctq is the 
Drude conductivity at a; = 0. 



Sjy = 4y4sem*crH£:^£:j, 
(41) 

^~Asem*aHE^E, (42) 



V. DISCUSSION AND CONCLUSION 

We have developed a semiclassical theory of non- 
linear transport and the photogalvanic effects in non- 
centrosymmetric media. It has been shown that the 
terms in semiclassical kinetic equations describing the 
electron motion associated with the Berry curvature and 
elementary shifts in the real space contribute to the lin- 
ear and circular photogalvanic effects as well as to the 
quadratic I-V characteristics (as a static limit of the 
LPGE). Previously Bloch et al.— studied the skew 
scattering mechanism of the LPGE. However, since it 
makes no contribution to the CPGE the Berry-curvature 
and shift related mechanisms certainly dominate the cur- 
rent sensitive to the chirality of the electromagnetic wave. 

It should be noted that we have considered linear-in- 
electric-field corrections to the collision integral due to 
the quantum coordinate shifts, see the term —eE{t) rpip 
in the (5-function of Eq. ([6|) . A comparable contribution 
to the tensors % and 7 comes from linear-in-£; corrections 

to the squared matrix element ^|V^'p|^^ allowed in non- 

centrosymmetric systems 1^ 

The following estimations summarize different contri- 
butions and allow one to compare their relative role in 
the formation of a dc current as a quadratic response to 
an ac electric field. The skew scattering contribution can 
be estimated as 

Jo 



el 



(43) 



where v and s are the characteristic electron velocity and 
energy (er 3> fi.), I = vt is the mean free path, and the 
dimensionless parameter of the scattering asymmetry is 
given by, see Eq. (PU)) . 



7-(cl,a) 
P 



(44) 



For static defects with a dipole moment d the asymmetry 
parameter is estimated^ as d/ (eas), where as is the Bohr 
radius lefi^ / {m* e^) , as is the static dielectric constant. 
If the dipole moment d is anomalously small one can 
take into account in the second term of the right-hand 
side of Eq. (|^D|) a combined scattering simultaneously 
involving a defect and an acoustic phonon. The result is 
formally obtained by the replacement of ■Ci™''^ in Eq. ((43|) 

^Csk^V(''ph£), where Tph is the electron momentum 
relaxation time due to scattering by acoustic phonons, 
?sk^^ is the dimensionless asymmetry parameter of the 
electron-phonon interaction including the deformation- 
potential and piezoelectric mechanisms, see, e.g., Ref . [26l. 

The Berry-curvature and shift-related contributions to 
LPGE, xb and Xsh, are given by the constants A and 77, 
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respectively, see Eqs. p4)) and (|33p . In semiconducting 
pyroelectric materials one can use a general estimation 
A 7] S,piroao/{m*Eg), where ao is lattice constant 
and ^piro is the dimensionless asymmetry parameter. As 
a result, we have 

Xb ~ Xsh CpiroRc(cri^)^. (45) 

For a wurzitc-type semiconductor one has from Eq. (j25p 
^piroOO ~ hQ/Ac . 

It follows from Eqs. (|43|) and ([45|) that, for compa- 
rable asymmetry parameters ■Ci™'''' ^ind ■^piro, the skew- 
scattering contribution to the LPGE prevails over the 
Berry-curvature and shift related contributions. 

On the other hand, the Boltzmann kinetic equation 
yields no contribution to the CPGE even if skew scatter- 
ing is included: 7sk = 0. The circular nonlinear current 
appears only if Fp and/or Vpip are nonzero in which case 



one has 

7b ~ Aem*lm{a^) , 7sh riem*lm{a^) . 

In the static limit both values 7b and 7sh as expected van- 
ish. At finite frequences, however, the Berry-curvature 
and shift-related effects provide the main contribution to 
the CPGE. 
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